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' Abstract. The single-pion production reactions pp div'^ , pp —> npn^ and pp ppyr" were measured at a 

QQ I beam momentum of 0.95 GeV/c {Tp ^ 400 MeV) using the short version of the COSY-TOF spectrometer. 

The central calorimeter provided particle identification, energy determination and neutron detection in 
addition to time-of-fiight and angle measurements from other detector parts. Thus all pion production 
channels were recorded with 1-4 overconstraints. The main emphasis is put on the presentation and discus- 
sion of the np-K^ channel, since the results on the other channels have already been published previously. 
The total and difi^erential cross sections obtained are compared to theoretical calculations. In contrast to 
the ppTr" channel we observe in the npvr"'" channel a strong infiuence of the A excitation. In particular the 
pion angular distribution exhibits a (3 cos^O + 1) dependence, typical for a pure s-channel A excitation 
and identical to that observed in the dix^ channel. Since the latter is understood by a s-channel resonance 
m the ^D2 pn partial wave, we discuss an analogous scenario for the pnn channel. 



■ PACS. 13.75.Cs - 14.20.Gk - 14.20.Pt - 25.10.+S - 25.40. Ep 



1 Introduction 

Single-pion production in the collision between two nu- 
cleons is thought to be the simplest inelastic process be- 
tween two baryons, nevertheless its understanding is still 
far from being satisfactory - both from the theoretical 
and the experimental point of view. In a recent publi- 
cation p] - in the following denoted by (I) - we have 
presented the first kinematically complete measurement 
for the pp ppK^ channel at a beam momentum of 0.95 
GeV/c (corresponding to Tp =397 MeV). Although A pro- 
duction in a relative s-wave is prohibited in this reaction 
channel, we have seen that the angular distributions indi- 
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cate the presence of significant A production in a relative 
p-wave already close to threshold. 

In this work we present our results for the pp — > pmr~^ 
channel at the same incident energy. We will show that 
here, where A production in relative s-wave is allowed, 
indeed this production process is overwhelmingly domi- 
nant and characterizes the differential observables in this 
channel. 

The reaction of interest here has been investigated in 
the near-threshold region already previously by a series of 
measurements at JINR p], KEK [SH], TRIUMF [5], CEL- 
SIUS |6] and notably at lUCF [7118119(110] . since there ex- 
periments have been carried out both with polarized beam 
and target. In these measurements very close to threshold 
it was already noted ^ that at beam energies of Tp = 
320 MeV, i.e. only 20 MeV above threshold, an onset of 
the A excitation is seen in the pion angular distribution. 



2 



S. Abd El-Samad et al.: Single-Pion Production in pp Collisions at 0.95 GeV/c (II) 




Fig. 1. Graphs used for the theoretical calculations. Top: t- 
channel approach (PWIA, direct part only), bottom: s-channel 
approach. 

The latter changes from isotropic at energies below 300 
MeV to anisotropic at 320 MeV suggesting there a 30 - 
40 % contribution of the A excitation. In the TRIUMF 
measurements, which were performed at beam energies of 
Tp = 420 and 500 MeV, it was noticed that at these en- 
ergies the pion angular distributions of pp dir'^ and 
pp — > prnr^ reactions are not only very similar, but also 
very close to the angular distribution expected from a pure 
A excitation. In terms of partial waves this means a pre- 
dominance of the pp partial wave at these energies 
as is in fact the outcome of the SAID [11] partial wave 
analysis for the pp dn^ reaction. Hence after presen- 
tation and discussion of our data we will confront them 
with theoretical i-channel calculations (Fig.l , top) as well 
as with s-channel calculations (Fig. 1, bottom). The lat- 
ter account for the striking partial wave dominance 
and assume for simplicity that all other partial waves give 
negligible contributions. In both cases an excitation of the 
A resonance is assumed. 



2 Experiment 
2.1 Detector setup 

Since the experimental setup has been discussed in detail 
already in (I), we give here only a short account. The mea- 
surements have been carried out at the Jiilich Cooler Syn- 
chrotron COSY using the time-of-flight spectrometer TOF 
at one of its external beam Hues. At the entrance of the de- 
tector system the beam - collimated to a diameter smaller 
than 2 mm - hits the LH2 target, which has a length of 
4 mm, a diameter of 6 mm and 0.9 fj,m thick hostaphan 
foils as entrance and exit windows. At a distance of 22 
mm downstream of the target the two layers of the start 
detector (each consisting of 1 mm thick scintillators cut 
into 12 wedge-shaped sectors) were placed followed by a 
two-plane fiber hodoscope (96 x 96 fibers, 2 mm thick 
each ) at a distance of 165 mm from target. Whereas the 
start detector mainly supplies the start time signals for the 
time-of-flight (TOF) measurements, the fiber hodoscope 
primarily provides a good angular resolution for the de- 
tected particle tracks. In its central part the TOF-stop 



detector system consists of the so-called Quirl, a 3-layer 
scintillator system 1081 mm downstream of the target - 
and in its peripheral part of the so-called Ring, also a 3- 
layer scintillator system built in a design analogous to the 
Quirl, however, with inner and outer radii of 560 and 1540 
mm, respectively. Finally behind the Quirl a calorimeter 
was installed for identification of charged particles and of 
neutrons as well as for measuring the energy of charged 
particles. The calorimeter consists of 84 hexagon-shaped 
scintillator blocks of length 450 mm, which suffices to stop 
deuterons, protons and pions of energies up to 400, 300 
and 160 MeV, respectively. The energy calibration of the 
calorimeter was performed by detecting cosmic muons. 

2.2 Particle identification and event reconstruction 

In the experiment the trigger required two hits in the Quirl 
and/or Ring associated with two hits in the start detec- 
tor. Tracks of charged particles are reconstructed from 
straight-line fits to the hit detector elements. They are 
accepted as good tracks, if they originate in the target 
and have a hit in each detector element the track passes 
through. In this way the angular resolution is better than 
1° both in azimuthal and in polar angles. If there is an 
isolated hit in the calorimeter with no associated hits in 
the preceding detector elements, then this hit qualifies as 
a neutron candidate (further criteria will be discussed be- 
low). In this case the angular resolution of the neutron 
track is given by the size of the hit calorimeter block, i.e. 
by 7 - 8°. By construction of the calorimeter a particle 
will hit one or more calorimeter blocks. The number of 
blocks hit by a particular particle is given by the track 
reconstruction. The total energy deposited by this parti- 
cle in the calorimeter is then just the (caHbrated) sum of 
energies deposited in all blocks belonging to the particular 
track. 

In order to have maximum angular coverage by the 
detector elements and to minimize the fraction of charged 
pions decaying in flight before reaching the stop detec- 
tors, the short version of the TOF spectrometer was used. 
In this way a total polar angle coverage of 3° < 6*'"^ < 
49° was achieved with the central calorimeter covering the 
region 3° < 0'"^ < 28°. For fast particles the energy res- 
olution of the calorimeter amounting to 4% is superior to 
that from time-of-flight measurements due to the short 
path length. However, the time-of-flight resolution is still 
much better than the AE resolution of the Quirl elements. 
Hence, for particle identification, instead of plotting AE 
versus Ecai, the uncorrected particle energy deposited in 
the calorimeter, we utilize the relation AE ~ {z/ (3)'^ with 
the particle charge z = 1 and plot 1 //3^ versus Ecai , where 
the particle velocity l3 — v/c is derived from the time-of- 
flight measurement. 

By identifying and reconstructing the two charged tracks 
of an event the exit channels dir'^ , npn^ and ppn'^ can be 
separated. Kinematically the maximum possible labora- 
tory (lab) polar angles are « 9° for deuterons and w 32° 
for protons (and neutrons). Hence 86% of the angular cov- 
erage for protons and neutrons from single pion produc- 
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tion are within the angular acceptance of the calorimeter. 
For charged pions the angular coverage has been much 
lower with this setup, since kinematically they can ex- 
tend up to 6*'"'' = 180°. Hence within the angular cover- 
age of Quirl and Ring the angular acceptance for 7r+ has 
been ^40% only. Nevertheless most of the phase space 
part necessary for a full coverage of the physics in sin- 
gle pion production has been covered (see below) by these 
measurements due to the circumstance that the center- 
of-mass (cm) angular distributions have to be symmetric 
about 90° because of identical collision partners in the 
incident channel. 



2.3 Selection of the npn^ channel 

The np7r+ channel is selected by identifying proton and 
pion in the calorimeter or only the proton in the calorime- 
ter, when the second charged track is in the Ring. In ad- 
dition, the missing pir mass MMp.^ has to meet the con- 
dition 900 MeV/c^ < MMp^ < 980 MeV/c^ . Also to 
suppress background from the d7r+ channel - in particu- 
lar when the deuteron breaks up and appears as a proton 
in the calorimeter - the pTr"*" track is required to be non- 
coplanar.To this end we determine the variable A<P, which 
is defined as the projection of the opening angle between 
two tracks onto the plane normal to the beam vector. That 
way we have always A<1> < 180°. The according histogram 
is displayed in Fig. 2. We see that p7r+ events stemming 
from deuteron breakup give rise to a large peak near 180°. 
In order to get rid of these events we introduce the con- 
straint Alp < 160°. From Fig. 2 we see that in prinicple 
a cut A<P < 170° would already be sufficient to eliminate 
this background. However, in order to be on the safe side 
and to have no perceptible tails from the pp dn'^ reac- 
tion in the final data sample, we used the more rigorous 
constraint. From Fig. 2 we see that this constraint does 
not cut away significant pieces of physics information 0. 
Any extrapolation into the cut region, be it by model or 
by phase space, will introduce uncertainties on the level 
of less than two percent. 

Aside from this physics background due to deuteron 
breakup we do not find any sizable background in the 
final data sample, as has been checked by control mea- 
surements, where the target cell was empty. 

Further on, the neutron 4-momentum is reconstructed 
from the 4-momenta of proton and pion and it is checked, 
whether a calorimeter block in the corresponding {0, 'P) 
region recorded a hit without any additional entries recorded 
in the preceding detector elements of the Quirl. If these 
conditions are met, a neutron track is assumed. That way 
On and <Pn are determined by the location of this calorime- 
ter block. By this method we obtain a neutron detection 
efficiency of 36 %. Thus having only the neutron energy 
undetermined experimentally we end up with 3 kinematic 
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^ we note that the data sample with the cut A$ < 170° 
leads to results, which are practically indistinguishable from 
those shown here 



Fig. 2. Distribution of the planarity angle A(j) defined in the 
center-of-mass system for the pp — + pmr^ reaction. Data of this 
work are shown by full circles and phase space by the shaded 
area. Solid and dashed lines denote s-channel and t-channel 
calculations, respectively, as discussed in the text. For ease 
of comparison the calculations have been normalized to the 
experimental total cross section. The constraint Afp < 160° 
for the elimination of events stemming from deuteron breakup 
is indicated by the vertical dash-dotted line. 



overconstraints for this channel. Corresponding kinematic 
fits were applied. 

The luminosity of the experiment was determined from 
the analysis oi pp elastic scattering, see (I). All data have 
been efficiency corrected by MC simulations of the de- 
tector setup by using the CERN GEANT3 [12] detector 
simulation package, which accounts both for electromag- 
netic and hadronic interactions of the ejectiles with the 
detector materials. 



2.4 Uncertainties 

The final data sample contains about 80000 good events, 
i.e. the statistical errors are on the percent level and hence 
of minor importance in comparison with systematical un- 
certainties. One source of systematic errors are beam align- 
ment and quality. The requirement that the beam has to 
pass a 2 mm collimator without significant halo assures 
good alignment and quality of the beam. The alignment 
is also verified by the fact that the angular distributions 
in the overall center-of-mass system have to be symmetric 
about 90 °, see next section. In addition polarized beam 
measurements with the identical detector setup used here 
also give no hint for noticeable misahgnments [13]. 

A much more severe source for systematic uncertain- 
ties concerns the acceptance and efficiency corrections of 
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the data. As discussed above the detector covers the com- 
plete momentum range , however, not the complete angu- 
lar range of the reaction of interest. This means, that the 
acceptance correction has to rely on extrapolations into 
unmeasured angular regions. From the Dalitz plots dis- 
played in Fig. 4 we see that - according to all we know 
from this reaction - these unmeasured regions are the ones 
with the lowest cross section. Hence systematic errors in- 
troduced by these extrapolations ought to be of minor im- 
portance. To quantify this statement we compare the data 
resulting from corrections with MC simulations, which are 
based either on pure phase space or alternatively on a 
model description of the pp — > npn^ reaction. A model, 
which is trimmed to describe all essential features of the 
data, is the appropriate tool for a reliable acceptance and 
efficiency correction by a MC simulation, which passes the 
ejectiles from the reaction of interest through the virtual 
detector. It has the potential of providing a selfconsistent 
procedure for these corrections. Contrary to this the pure 
phase space description of the reaction is - as we easily 
can see from the experimental results displayed in Figs. 2 - 
6 - inadequate though convenient. It may serve, however, 
as a very conservative estimate of the systematic uncer- 
tainties due to the acceptance and efficiency corrections. 
As examples we display in Fig. 3 our data for the invari- 
ant pn mass Mp„ and for the pion angular distribution in 
the center-of-mass system (j{0^'^^) evaluated (properly) 
by use of a model, which fits the final data (solid circles), 
and alternatively by use of pure phase space (open circles), 
respectively. The use of the latter for the extrapolation 
into unmeasured regions brings the data, of course, some- 
what closer to the phase space predictions. However, de- 
spite the fact that phase space and model predictions are 
vastly different, the effect of the inadequate phase space 
correction on the data is still very moderate. A reaHstic 
upper limit for systematic uncertainties due to acceptance 
and efficiency correction will be much smaller than the 
differences between open and soHd circles in Fig. 3. As a 
realistic estimate for this kind of systematic uncertainties 
we conclude that they are within the size of the symbols, 
which are used in Figs. 2 - 6 for displaying our experi- 
mental results. We note that the use of either t-channel 
or s-channel calculations in the MC simulations does not 
lead to any noticeable differences in the acceptance and 
efficiency corrected data. 

3 Results 

Due to the identity of the coUision partners in the entrance 
channel the angular distributions in the overall center-of- 
mass system have to be symmetric about 90 °, i.e. the full 
information about the reaction channels is contained al- 
ready in the interval 0° < 6*'^™ < 90° . Deviations from this 
symmetry in the data indicate systematic uncertainties in 
the measurements. Hence we plot - where appropriate - 
the full angular range, in order to show the absence of 
major systematic errors in our measurement. 

The total cross section of 0.47(2) mb for the pp 
pn-K^ reaction at Tp = 400MeT^ has already been given 
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Fig. 3. Differential cross sections in dependence on the in- 
variant mass Mpn (left) and the pion scattering angle O'^^' 
(right) in the center-of-mass system for the pp prnr^ re- 
action. The full circles show our data properly corrected by 
a self-consistent MC simulation, i.e. using a reaction model, 
which is in agreement with the final data. The open circles 
derive from MC simulations using just pure phase space. The 
phase space distributions for the differential spectra are shown 
by the shaded area. Solid and dashed lines in the right figure 
show Legendre fits to the filled and open circles yielding a2 = 
0.96(2) and =0.79(2), respectively. 



in (I) . It is roughly a factor of two smaller than that for 
the pp (i7r+ channel and five times larger than that 
for the pp ppn^ channel, see Table 1 in (I). For the 
comparison with data at neighbouring energies [T4ll3l [5ll8] 
see discussion in the next section. 

Differential distributions are shown in Figs. 3-6. Since 
a 3-body system in the final state has five independent 
variables, we choose to present the three invariant mass 
spectra Mp„, Mpjr+ and M„7r+ as well as the proton and 
pion angular distributions in the overall center-of-mass 
system (cms). Because of the Hmited acceptance for pi- 
ous only the angular range 0° < 6*^™ < 90° is covered. 
However, as pointed out above, due to the required 90° 
symmetry of the angular distributions the full informa- 
tion is contained in the data. 

In contrast to the ppn'^ channel the Dalitz plots for 
the pmr^ channel are far from being fiat - see Fig. 4. This 
is also borne out in the projections, the invariant mass 
spectra Mp^+ and M„jr+ (see Fig. 5, top), which peak at 
the highest kinematically available masses. In contrast to 
these the Mp„ spectrum (Fig. 5, bottom) peaks at the 
lowest masses. Note that the high-mass enhancement is 
similar in both Mp^+ and M„^+ spectra. This is in agree- 
ment with the trend observed in bubble chamber data at 
Tp = 432 MeV 0. 

Whereas the proton angular distribution is close to 
fiat, the pion angular distribution is strongly anisotropic. 
It is fitted very well by a pure (3cos^6>^™ -I- 1) distribution 
(solid line in Fig. 3, right). 

Our results are in good agreement with the results 
from PROMICE/WASA j6], which are the only other un- 
polarized data available at Tp w 400 MeV. Note, however, 
that those measurements were carried out in a very limited 
phase space region only. E.g., the pion angular distribu- 
tion was only measured for cos{0'^"^) > 0.79. 
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Fig. 4. Dalitz plots for the invariant mass combinations 



versus M„ 



(left) and M^^+ versus M^^+ (right) as obtained 
for the pp pmr~^ reaction: data are shown on top and the MC 
s-channel calculation (see text) at the bottom. Note that the 
plots for the data are efficiency but not acceptance corrected. 
The tiny deviations from the elliptic circumference at the upper 
corners are due to the excluded beam-hole region. 
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Fig. 6. Angular distributions of protons (left) and pions (right) 
in the center-of-mass system for the pp — > pnn^ reaction. Data 
of this work are shown by full circles, WASA/PROMICE re- 
sults [6] by open circles (renormalized to atot of this work) and 
phase space by the shaded area. Solid and dashed lines denote 
s-channel and t-channel calculations, respectively, as discussed 
in the text. For ease of comparison the calculations have been 
normalized to the experimental total cross section. Given are 
also the Legendre coefficients a2 from a Legendre fit according 
to eq. 2 in (I) tgether with statistical errors. 
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Fig. 7. Differential cross sections in dependence on invariant 
masses Mpp (left) and M^^o (right) for the pp ppn^' reaction 
(from (I)). Phase space is shown by the shaded area. For the 
explanation of the curves see (I). 



4 Discussion of Results 
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Fig. 5. 

masses 



Differential cross sections in dependence on invariant 



and M„^+ for the pp — > pnTr""" reaction. 



Data of this work are shown by full circles and phase space by 
the shaded area. Solid and dashed lines denote s-channel and 
t-channel calculations, respectively, as discussed in the text. 
For ease of comparison the calculations have been normalized 
to the experimental total cross section. 



4.1 Comparison with the ppir^ channel 

For the comparison of prnr^ and ppn'^ channels we show 
once more results in Figs. 7 and 8, which we obtained 
previously for the ppir'^ channel, see (I). 

The proton angular distributions both in pri7r+ and 
pp7r° channels are close to isotropic (see Figs. 6 and 8, left). 
They nevertheless show a tendency of opposite curvature, 
which just may be some reflection of the very different 
pion angular distributions in both channels. At these low 
energies the pion angular distributions are described con- 
veniently by the Legendre ansatz 

ct{9^J) ~ 1 + 02 * {Scos^e^J - l)/2, 

where the parameter 02 is fltted to the data. This 
ansatz has also been used in (I) for the analysis of the 
pp7r° channel - see eq. 2 in (I) . Other notations related to 
our a2 parameter are as follows: in the notation of Ref. 
[5], Table III we have 02 = A2/AQ, in the notation of Ref. 
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Fig. 8. Angular distributions of protons (left) andpions (right) 
in the center-of-mass system for the pp ppn'^ reaction (from 
(I)). The solid curves represent Legendre fits, see (I). 



[H], Table V we have 02 = B/A and in the notation of Ref. 
|10) , Table II we have 02 = 2 * 600 (since ago = 1 there) . 

For the pion angular distribution in case of ppn'^ we 
find a slightly concave shape (Fig. 8, right; 02 = -0.12), 
whereas we get a strongly convex shape in case of the 
pmr'^ channel (Fig. 6, right; 02 = +0.96(2)). 

A striking difference to the ppir'^ channel is the dom- 
inance of the A excitation in the pmr~^ channel. In the 
ppir'^ channel A production in a s-wave relative to the ac- 
companying nucleon is prohibited by quantum numbers. 
However, as demonstrated in (I), the A production in rela- 
tive p-wave is allowed and indeed contributes significantly 
to this channel already close to threshold. 

The signature of A production in the pmT~^ channel 
is seen in all differential observables shown in Figs. 3-6. 
The Mpn spectrum (Fig. 5, bottom) is markedly different 
from the Mpp spectrum in the ppn^ channel (Fig. 7, left). 
Though the pn final state interaction (FSI) is somewhat 
larger than the pp FSI, it cannot account for the huge 
difference between both spectra at low masses. It rather 
is a refiection from the complementary spectra Mp^+ and 
M„^+, which peak just at the highest available masses 
(Fig. 5, top) in contrast to the more phase-space like Mp^o 
spectrum observed in the ppn'^ channel (Fig. 7, right). The 
high-mass enhancement in turn is most likely associated 
with A excitation. 



4.2 Pion angular distribution 

The most striking feature in our pmr~^ data at Tp = 400 
MeV is the (3 + dependence of the pion angular 

distribution. It is exactly this dependence, which is seen 
in ttN scattering in the A resonance region and which is 
also observed in the pp d7r+ reaction over a large range 
of incident energies. We note that also the pion angular 
distributions at Tp = 420 and 500 MeV, deduced from 
inclusive measurements of the pp — > prn:^ reaction in Ref. 
[5] are basically in accord with this dependence. 

We also note in passing, that in contrast to the sit- 
uation in the ppn^ channel, where we observe a strong 
dependence of the pion angular distribution on the rel- 
ative momentum q between the two nucleons in the exit 



channel, we do not find such a behavior here. The pion an- 
gular distributions stay essentially the same for different 
q-regions. This again is in support of a single dominant 
mechanism proceeding via a single partial wave. 

Whereas we find good agreement in the anisotropy of 
the pion angular distribution with the ones measured at 
TRIUMF at Tp = 420 and 500 MeV 0, we face a pro- 
found discrepancy with lUCF results There kinemat- 
ically complete measurements with both polarized beam 
and target are presented for Tp = 325, 350, 375 and 400 
MeV together with Legendre polynomial fits, which repro- 
duce the measured angular distributions of the polariza- 
tion observables. In Table II of Ref. [lO] the coefficients 
resulting from these fits are listed. This list also contains 
the anisotropy parameter 600 = 02/2 for the unpolarized 
differential cross section a{0^^). Since that paper con- 
centrates on the polarization observables, it does not di- 
cuss the unpolarized (t{0^'^'^) data in any detail. It is just 
noted that these have been deduced from p + n coinci- 
dences. As a result of their Legendre analysis they obtain 
a2 = 0.39 (9) for Tp = 400 MeV, a value, which is less than 
half of our value. The difference amounts to six standard 
deviations. We face here perhaps a similar problem as dis- 
cussed in (I) for the situation of the experimental data 
with the reaction pp ppn^ at the same beam energy 
of Tp = 400 MeV. There also a large discrepancy exists 
between the asymmetry parameter 02 found in the lUCF 
analysis and the one we found at COSY-TOF. Our result 
has meanwhile been confirmed by CELSIUS/ WAS A data 
[15] . In both these works it has been shown that previous 
results - like the one obtained at PROMICE/WASA [16] 
- suffered from insufficient phase space coverage combined 
with inappropriate extrapolations into the unmeasured re- 
gions. 

We note the following surprising trend in the results of 
the lUCF analyses concerning the 02 parameter. In agree- 
ment with previous lUCF work very close to threshold [8], 
where a strong increase in the a2 value with energy is ob- 
served (02 = 0, 0.12 and 0.31 at Tp = 294, 299 and 319 
MeV, respectively), a value of 02 = 0.34(7) at Tp = 325 
MeV is reported in Ref. [T^, Table II. This increase just 
refiects the increasing importance of the A excitation in 
the pp — > npK~^ reaction. However, above Tp = 325 MeV 
Ref. [I0j(see again Table II) finds that the 02 parameter 
no longer increases, but rather saturates at a value of 0.39. 
This implies that the A contribution also suddenly satu- 
rates despite the fact the A pole is still far from being 
reached at these energies. Also the .liilich model calcula- 
tions [17J, which are used in Ref. [10] for comparison with 
the polarization data predict a strongly increasing impor- 
tance of the A excitation at these beam energies - as also 
expected intuitively. 

In Ref. [1^ the discrepancy between their 400 MeV re- 
sult and the TRIUMF data at Tp = 420 MeV [5] was not 
discussed. The TRIUMF measurements were performed 
as a single-arm experiment with a magnetic spectrometer, 
which measures a 3-body reaction only inclusively. Such 
measurements are known to be very reliable for the de- 
termination of single-particle angular distributions, since 
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exactly the same detector is used for the particle detection 
at all angles minimizing thus the problem of acceptance 
and efficiency corrections. 

To summarize the discussion about 02 we face a severe 
discrepancy between the lUCF and COSY-TOF results at 
Tp = 400 MeV for both the pp — > ppir'^ reaction and the 
pp np-K^ reaction. However, there is good agreement of 
our result for the first reaction with the corresponding one 
from CELSIUS/WASA and for the second reaction with 
the one obtained at TRIUMF. We finally note that the 
uncertainty of 02 quoted in Fig. 6 only includes the statis- 
tical uncertainties of the data. If we include an estimate 
on systematic uncertainties (as discussed above in the sec- 
tion 2.4), we end up with an estimated total uncertainty 
of Aa2 « 0.08. 

The close similarity in the angular distributions of di:^ 
and pniT^ channels is in agreement with the predictions 
of Faldt and Wilkin p[8|, according to which the meson 
angular distributions coincide closely in bound state and 
breakup channels near threshold. 

The (3 -\- 1) dependence in the pp dw^ re- 

action is due to the ^1)2 partial wave in the incident pp 
channel combined with the constraint of a p-wave pion rel- 
ative to the deuteron in the exit channel. As known from 
phase shift analyses of this reaction [TT|ll9j this ^D2P par- 
tial wave dominates the pp dn'^ reaction basically from 
threshold up to Tp « 900 MeV. It is responsible for the 
resonance-like energy dependence of the total cross sec- 
tion and performs a perfect looping in the Argand dia- 
gram [11]. Thus this partial wave possesses all features to 
qualify for a s-channel resonance with I{J^) — 1(2'*') in 
the incident pp channel, the final c?7r+ channel and the in- 
termediate NA channel - see also discussion in Ref. [27j. 
From the total cross section of the ^D2P partial wave in 
the pp — > diT^ reaction [11] we see that the resonance en- 
ergy is some 60 MeV below the nominal NA threshold, 
whereas the resonance width of « 110 MeV corresponds 
to that of the A resonance. 

The experimental observation that up to Tp = 500 MeV 
essentially a (3 cos^0'^+ +1) dependence is observed also 
in the pp pmr^ reaction points to the predominance 
of the ^D2P partial wave in this reaction, too, as was 
already noted in Ref. [5]. From this also follows that the 
pn system at these energies is preferably in the deuteron- 
like quantum state I{J^) — 0(1+). 



4.3 Total cross section 

The energy dependence of the total cross sections for the 
pp dn~^ and pp — > np7r+ channels is displayed in Fig. 9. 
For the pp prnr^ channel at Tp = 400 MeV there are no 
experimental data to compare with. However, our result 
fits well to the trend given by the experimental results 
|14[[3ll5|[8] at neighbouring energies. In order to discuss in 
the next section the energy dependence in some broader 
context with regard to the dominance of the ^D2 partial 
wave, we plot in Fig. 9 the energy range from threshold up 
to Tp = 1.5 GeV, i.e. over the full region of A excitation. 



- 


s+t 


- 






























ir ■ \ ,5 □ 












d ; ■ N S 







0.2 0.4 0.6 0.8 1 1.2 1.4 



Tp [GeV] 



Fig. 9. Energy dependence of total cross sections for the reac- 
tions pp dTT"*" (squares) and pp pnn^ (circles) . The filled 
symbols denote results of this work, open symbols are from [HI 
[SUl]. Dash-dotted and dashed lines represent s-channel calcu- 
lations for d7r+ and pnTr"*" channels, respectively, normalized to 
the data at Tp = 400 MeV. The ^-channel calculation for the 
pmx^ channel is given by the dotted line and normalized to the 
data at Tp > 1 GeV, where it is known to provide a reasonable 
description. The solid curve denoted by s + t is the incoherent 
sum of both processes for the pp pmr'^ channel. 

4.4 Comparison to t- and s-channel calculations 

Finally we confront the data with simple theoretical calcu- 
lations according to the graphs depicted in Fig. 1. In this 
experimental paper it is not our aim to compare our data 
with complex calculations, which are beyond the scope of 
this work. We rather want to use these calculations in or- 
der extract the main physics message residing in the data. 

Since we have seen in the discussion above that A exci- 
tation and decay is the dominant reaction mechanism and 
^I?2 is by far the dominant partial wave at least in the 
energy region Tp = 400 - 500 MeV, we expHcitly neglect 
any non-resonant terms in the reaction process. 

We start with a (properly antisymmetrized) i-channel 
approach in first-order impulse approximation, which is 
known to provide reasonable results for the angular dis- 
tributions at energies Tp > 1 GeV [20ll2Tll22ll23l[24] . For 
our energy we show such a t-channel calculation by the 
dashed lines in Figs. 2, 3, 5 and 6. We see that the pion 
angular distribution is not reproduced correctly. The rea- 
son for this failure is simple. Since in this approach the A is 
excited by pion exchange between the two colliding nucle- 
ons, the reference axis for the pion decay of the excited A 
resonance is the momentum transfer q and not the beam 
axis. Note that the latter, however, is the reference axis 
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Fig. 10. Energy dependence of total cross sections as calcu- 
lated with our s-channel approach for the pp — > dir^ (dash- 
dotted) and pp —>■ pn-K^ (dashed cureve) reactions, normalized 
to the data points from this experiment at Tp = 400 MeV (filled 
square and circle, respectively; see also Fig. 9). For comparison 
the short dash-dotted and dotted lines give the SAID [11] so- 
lutions SP07 and SP40, respectively, for the partial wave 
in the npn^ channel (see text) , whereas the solid curve shows 
the SAID solution SP93 for the partial wave in the d-K'^ 
channel. 



for angles in the center-of-mass frame. Since the momen- 
tum transfer varies in dependence of the nucleon scatter- 
ing angles, which are integrated over in the pion angular 
distribution, the calculated intrinsic (3 cos^O'^^ + 1) de- 
pendence in the ^-channel approach appears to be smeared 
out in the cms pion angular distribution. 

In order to establish the beam direction as the appro- 
priate reference axis for the pion angular distribution in 
the calculation one either has to reiterate the pion ex- 
changes and sum them up to infinity or simply make a 
s-channel ansatz. The latter is easily made by a Breit- 
Wigner ansatz for the resonance (Fig. 1, bottom), 
which dissociates into a NA system in relative s-wave fol- 
lowed by the decay of the A resonance. This grants the 
observed (3 cos^©^™ -1-1) dependence for the pions as well 
as the desired shape of the invariant mass spectra. To sim- 
ulate the pp dw^ reaction within this ansatz we impose 
the Hulthen wave function as condition for the Fermi mo- 
mentum between proton and neutron. For the pp —>■ pm:^ 
reaction we impose a FSI interaction of Migdal-Watson 
type |25[[26] on the pn system. These calculations (nor- 
malized to our results for the total cross sections at Tp = 
400 MeV) are shown by the Dalitz plots in Fig. 4, bottom 
and by the solid lines in Figs. 2,3, 5 and 6. They give a 



good account of the differential distributions both for the 
invariant masses and for the particle emission angles. 

Mass (2.07 GeV) and width (160 MeV) of the /(J^) = 
1(2+) resonance in this s-channel ansatz have been chosen 
to reproduce the energy dependence of the total cross sec- 
tion of the pp rfTT"*" reaction in the region of the ^D2P 
dominance, i.e. up to Tp « 700 MeV. That way we simulta- 
neously get a good description for the energy dependence 
of the total cross section for the pp pmr'^ reaction up to 
Tp w 600 MeV. If combined with the amplitude of the t- 
channel approach we even obtain a reasonable description 
of the total cross sections up to the GeV region (Fig. 9). 

Finally we compare in Fig. 10 our simple s-channel 
calculations for the ^I?2 partial wave in the np7r+ channel 
to the results of SAID [H] partial wave analyses of NN 
scattering and pp dTr"*" reaction. The imaginary part of 
the partial wave, which is obtained from the analysis 
of elastic NN scattering, gives rise to the total inelastic 
cross section (7ineiCT)2) for the ^D2 partial wave, which 
in turn is given by the sum of the pion-production cross 
sections in this partial wave. Since in the pp ppi:^ reac- 
tion the ^D2 partial wave is highly suppressed due to the 
suppression of the A excitation, we have in good approx- 
imation 

^pp^npTT+i D2) ~ (Jineli D2) ~ Crpp^d7r+ { D2) 

The thus derived results are shown in Fig. 10 for the 
SAID solutions SP07 and SP40 from NN scattering anal- 
yses together with the SAID solution SP93 for the ^D2 
partial wave in the fi7r+ channel. We see that our s-channel 
ansatz overpredicts the ^D2 part of the pp — > dTr"*" cross 
section at higher energies as expected, since we have as- 
sumed for simplicity that only this partial wave contributes. 
For the np7r+ channel the SP07 solution underpredicts 
the cross section in the near-threshold region strongly, 
whereas the SP40 solution, which is made for the region 
Tp < 400 MeV, provides a proper energy dependence, how- 
ever, is in the absolute scale substantially above the data 
and the s-channel calculation in this energy region. At 
higher energies up to 800 MeV our calculation is essen- 
tially between both solutions. This comparison with the 
SAID solutions for the ^D2 partial wave shows that our 
simple model description is not in severe contradiction to 
the SAID results, where the large spread between both 
solutions points to still substantial ambiguities in the NN 
partial wave analyses concerning the imaginary parts of 
partial waves. 



5 Summary 



priTT'^ re- 



We have presented measurements of the pp 
action at Tp w 400 MeV. In this energy region they are 
the first exclusive ones of solid statistics to cover most 
of the reaction phase space. The differential distributions 
are characterized by the overwhelming dominance of the 
A excitation. The pion angular distribution coincides with 
that observed in the pp (i7r+ channel, which in turn is 
dominated by the resonating ^I?2-P partial wave. From 
this coincidence we conclude that also the pn7r+ channel 
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is dominated by the same s-channel resonance in the near- 
threshold region. The correlation between boundstate and 
breakup channels is in accordance with the predictions by 
Faldt and Wilkin ^18]. 
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discussions with R. A. Arndt, L. Alvarez-Ruso, D. Bugg, C. 
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